% Расчет ДН параболической антенны
Angls=720;
delta=2*pi/(12*Angls);
teta_p=-pi/12:delta:pi/12-delta; % Полярный угол параболической антенны
F_parab=0;
for kk=1:Angls/2
  U(kk)=(k*R0*sin(teta_p(kk)));
  L1=2/(U(kk))*besselj(1,abs(U(kk)));
  L2=8/(U(kk))^2*besselj(2,abs(U(kk)));
  L3=48/(U(kk))^3*besselj(3,abs(U(kk)));
  F_parab(kk)=((1+a2+a4)*L1-(a2/2+a4)*L2+(a4/3)*L3);
end
% Нормировка матрицы 
maxF_parab=max(abs(F_parab));
F_parab=abs(F_parab./maxF_parab);
% Симетричная половина ДН
for ij=1:Angls/2-1
  F_parab(Angls/2+ij)=F_parab(Angls/2-ij);
end
for ij=1:Angls-1
  F__parab(ij+1)=F_parab(ij);
end
F_parab=F__parab;
clg;figure(gcf); 
plot( teta_p(1:(Angls-1))*180/pi,(F_parab(1:(Angls-1))),'m-');
title('Диаграмма зеркала');
xlabel('полярный угол ,в грусах');
ylabel('величина нормированной ДН');
grid on
pause
visual(teta_p(1:Angls-1), abs(F_parab(1:Angls-1)), Angls-1); % множитель системы
  title('ДН параболической антенны в полярных коорд.')
xlabel('полярный угол ,в грусах');
ylabel('величина нормированной ДН');
pause
close
 
